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In the present work, utilizing the Crowell and 
Labuda model, the MTF degradation due to lateral 
diffusion is calculated for a back-illuminated CCD 
imager for typical device parameters. The discrete 
nature of the CCD and finite size of the photosensi- 
tive elements, with its less than ideal collection effi- 
ciency for photogenerated charge, result in an additional 
MTF degradation. The Fourier transform approach 
is utilized to calculate the effective point-spread func- 
tion for these processes in the time domain. Experi- 
mental data is presented on the point- spread function 
for a three-phase, double-level anodized aluminum 
160 X 100 thinned and back-illuminated CCD imager 
and compared with the theoretical results. A simple 
modification of the Crowell and Labuda model sug- 
gested by these results is presented. Finally, the 
effects of the less than ideal charge transfer efficiency 
in practical devices are calculated and applied to pre- 
dict the anticipated MTF and point- spread function for 
a 400 X 400 large-area CCD imager currently under 
development. 

I. INTRODUCTION 

The charge -coupled device offers a new and significant approach to 
achieving the once elusive goal of a viable large-area solid-state imager 

This paper presents the results of one phase of research carried out at the Jet 
Propulsion Laboratory, California Institute of Technology, under Contract No. 
NAS 7-100, sponsored by the National Aeronautics and Space Administration. 
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possessing the format size and element density comparable to a conventional 
vidicon. Considerable progress has been made towards achieving this goal, as 
reflected by the performance data presented at this Symposium on large-area 
array imagers. This paper discusses the MTF, point-spread function, and the 
effects of charge transfer inefficiency for such large-area arrays. Expressions 
for the MTF and point-spread function are derived and applied to a thinned, 
back-illuminated 400 X 400 area imager (Ref. 1). The results indicate the 
anticipated performance levels for such a CCD and the significant role of 
charge transfer efficiency for large-area arrays* 

II. MTF AND POINT-SPREAD FUNCTION FOR LATERAL DIFFUSION 

Figure 1 illustrates the simplified model utilized by Crowell and Labuda 
(Ref. 2) to calculate the quantum efficiency r lg and the MTF for a silicon diode 
array target. The depletion region of width L^ - L a is assumed to be totally 
depleted and a region of zero lateral conductivity, tIq and the MTF are calcu- 
lated by evaluating, for a sinusoidal optical input, the steady-state diffusion 
current 


J = 
P 



( 1 ) 


flowing into the depletion region from the field-free region and adding the photo- 
generated current from the direct photon absorption in the depletion region. J 

r 

is evaluated from the solution of the continuity equation 


-DV 2 p + “F = G(x, y, z) 


( 2 ) 


subject to the boundary conditions 


= If 


y = 0 


(input surface) 


P = 



(depletion region edge) 


(3) 
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where 


p = minority carrier density in excess of thermal equilibrium 
t = minority carrier lifetime 

2 

D = minority carrier diffusion constant (cm / sec) 

G(x, y, z) = minority carrier generation rate per unit volume 
S = surface recombination velocity (cm/sec) 

The photogeneration rate G(x, y, z) is assumed to vary sinusoidally in the trans- 
verse direction and exponentially in the y direction, and is given by 


0 

G(x, y, z) = -^-(1 - R)a exp (-ay) (1 + cos kx) (4) 

where 

Nq = peak incident photon flux 

= silicon absorption coefficient at the wavelength of interest 
R = reflectivity of silicon at the wavelength of interest 
k = Ztt/X = spatial frequency of the optical input 

The resulting solution can be written in the form 


N 0 

Jp(x) = — (ti 0 + n k COS kx) 


with 


= »L(1 - R) f2(»L -f SL/D) - - gj exp (-«L„) _ _i_ , .1 

a 2 ! 2 - 1 L P + + K « L < a a j 


-(1-R) exp (-ffL b ) 


(5) 
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and 


■where 


n o = 


k=0 


P* = ( 1± 'lB ex I> ± ir 

Lq = diffusion length = (Dt)^^ 


L = 


(l + k 2 L 2 ) 1/2 


The MTF, defined to be the peak-to-peak minority carrier current entering the 
continuous p region normalized with respect to the response for uniform light, 
is simply the ratio "H /"H q • The resulting MTF can be written as a sum of a 
frequency-dependent component and a constant, 


_k _ r|(a>) 



(6) 


where 


C = exp (-<*L a ) {l - exp [-<*(1^ - L a )}} 


r l( w )/ r lg is the output modulation due to absorption in the field-free region and 
decreases due to lateral diffusion as the input spatial frequency is increased. 
C/tIq, on the other hand, is equal to the fraction of the total collected charge 
contributed by direct photon absorption in the depletion region and contributes 
a constant output modulation. This latter term arises from the total depletion 
approximation. 

The above results can also be applied to a thinned, back-illuminated CCD 
imager such as that shown schematically in Figure 2. For a typical oxide thick- 
ness, clock voltages, and substrate resistivity, the maximum depletion depth is 
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approximately 5 pm. The depletion depth x^, however, is not constant throughout 

the unit cell but varies schematically, as illustrated in Figure 2. Typically, for 

a three-phase CCD, only approximately 1/3 of the unit cell is depleted to the 

maximum depth, with the remainder of the unit cell being only partially depleted. 

The exact profile can be obtained from the surface potential profile <p (x), since 

s 



where e is the dielectric constant of silicon and N_ is the acceptor density of the 
substrate. Additionally, (p (x) is, in general, anisotropic, with the doping and 
width of the channel stops being important factors in the direction perpendicular 
to the channel stops. However, the quantum efficiency is relatively insensitive 
to the detailed shape of the potential wells. Although partial depletion introduces 
a transition region of charge collection uncertainty and cross modulation between 
the unit cells, the total collection efficiency remains unity since the recombina- 
tion loss is small in this region. 

Finally, at wavelengths greater than 0, 8 pm, one can expect interference 
effects due to multiple reflections from the structured side of the CCD. Seib 
(Ref. 3) has extended the Crowell and Labuda model to include the effects of 
multiple reflections for normal incidence. His results show that for wavelengths 
greater than 0. 8 pm, although the calculated 9 q is increased over the single pass 
case, the resulting MTFs are nearly identical. Thus, although the total deple- 
tion condition is not realized in practice by a back-illuminated CCD, and multiple 
reflections are not considered at the longer wavelengths, one nevertheless should 
obtain a reasonable estimate of the quantum efficiency t^q and the MTF degrada- 
tion due to lateral diffusion from the single-pass Crowell and Labuda model. 
Equation (6) has been evaluated utilizing typical CCD parameters for two sub- 
strate thicknesses, and the results are shown in Figures 3 and 4. For a thinned 
array, the results are relatively insensitive to the value of Lg and to S. The 
key parameter, especially at short wavelengths, is the thickness of the field- 
free region. MTF (®) is the constant to which the curves converge and is equal 
to C/t1q. 


196 



The results can also be expressed in the time domain by calculating the 
inverse Fourier transform of Equation (6) to yield the point-spread function for 
lateral diffusion. 


PSF(x) = <^“ _1 1^1 = exp (-iwx)du + 6(x) (8) 

Since the lateral diffusion process possesses even symmetry, it is only neces- 
sary to evaluate the real part of the complex Fourier transform. Equation (8) 
was evaluated by numerical integration. The results are shown in Figures 5 
and 6 and represent the spatial charge spread (impulse response) of the field- 
free region to a delta optical input. 

HI, EFFECTS OF SAMPLING AND THE FINITE PIXEL SIZE 

In addition to the MTF degradation resulting from lateral diffusion, the 
discrete and finite aperture of a CCD imager imposes a further geometric limi- 
tation on the MTF. The resulting MTF loss can be estimated by assuming total 
depletion and dividing the depletion region into contiguous pixels of width d, 
where d is the pitch of the pixels, as shown in Figure 2. Under these simplify- 
ing assumptions, the MTF for the sampling and collection process is given by 


Rp * sine y (9) 

The total pixel MTF becomes r lj c / r lQ x Rp and is plotted in Figure 7 for a typical 
pixel size of 1 mil. Note that in the spatial frequency range of interest (f < ^yq^' 
the pixel size is the dominant factor. Additionally, the spectral dispersion in 
the MTF curves has been reduced. The corresponding point-spread function has 
been calculated by first truncating the MTF expression to zero at f = 4000 c/mm 
and numerically integrating up to this limit. The results are shown in Figure 8. 
The error introduced by truncation is negligible for \ = 0. 4 and 0. 6 pm, 
whereas for X = 0.8 pm, the actual point- spread function should be slightly 
more rectangular than calculated. 
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The response of a single pixel is given by 


R(x) = K I I(x')PSF(x» - x)dx» 


( 10 ) 


where I(x') is the spread function for the optical input. If I(x') is narrow in com- 
parison with PSF(x), it follows from the sifting properties of the delta function 
that 


R(x) S PSF(x) 


(ID 


Thus, by measuring the response of a single pixel as a point or narrow slit light 
source is scanned across the pixel, the point-spread function of the device can 
be experimentally determined. The next section will describe the results of 
such a measurement on a thinned, back-illuminated CCD imager. 

IV. EXPERIMENTAL RESULTS 

The characteristics of the tested CCD and the operating conditions for these 
measurements, respectively, are summarized below. Design and general per- 
formance data has been presented previously on this type of CCD imager (Ref. 4). 

(1) CCD imager characteristics 

(a) Three-phase, double-level anodized aluminum 

(b) Surface channel 

(c) Thinned and back-illuminated 

(d) Array size - 160 pixels horizontal by 100 pixels vertical 
(160 x 120 mils) 

(e) Pixel size - 1.0 X 1.2 mils 

(2) Operating conditions for pixel sensitivity profile measurements 

(a) Serial clock frequency = 100 kHz 

(b) Serial fat zero » 10% 
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(c) Exposure time = 250 msec (shuttered) 

(d) Total frame time = 423 msec 

The point-spread function measurements were made with the apparatus 
shown in Figure 9. The point light source was a 5 -mil Xenon arc lamp, which 
was imaged through a 10 x F/2 microscope objective. The CCD was approxi- 
mately 200 cm from the light source. Spectral and neutral density filters were 
inserted into the optical path to vary the spectral pass band and intensity of the 
light spot. The light spot could be viewed simultaneously through a microscope 
eyepiece via a beam splitter as it was being focused on the CCD input surface. 
The diameter of the minified light spot, as visually determined by focusing the 
light spot onto a stage micrometer placed in the same plane as the CCD input 
surface, was less than 4 pm. The light spot was stepped across the CCD in the 
parallel transfer section in 2. 2-pm increments in the image plane. The output 
from a 20 x 20 area (400 pixels) was recorded after a 12-bit A/D conversion by 
a digital tape recorder after each 2. 2-pm increment. The data was taken from 
a computer printout of the 20 x 20 area. Figure 10 shows the typical video out- 
put as the spot was scanned across three adjacent pixels. The assymetric 
response is due to the redistribution of a fraction of the charge from the primary 
pixel into the trailing pixels because of charge transfer inefficiency. 

Typical results showing the simultaneous output from three adjacent pixels 
as a function of spot position are shown in Figures 11 and 12, respectively. A 
Beck glass filter with a long-wavelength cutoff of 0. 45 pm was utilized. The 
spot intensity was adjusted to be approximately one-half of full well. The spot 
was positioned as close as possible to the output to minimize charge transfer 
efficiency effects. The peak amplitude increased from pixel to pixel because of 
the localized shading in the region near the output. The slight assymetry in the 
measured profiles is due to charge transfer inefficiency. The solid lines enclos- 
ing the sensitivity profiles are the boundaries of the respective pixel. Similar 
measurements were also made with a Beck red glass filter, and essentially 
identical results were obtained. 

Comparison of the sensitivity profiles with the calculated point-spread 
functions of Figure 8 shows that the measured profile is narrower at the center 
and exhibits more response outside the pixel than predicted by the Crowell and 


199 



Labuda model, indicating some pixel overlap. A simple modification of this 
model is suggested by these results, as shown in Figure 13. The collection 
regions defining the pixels are assumed to be overlapping and trapezoidal in 
shape to account for the less than ideal isolation between pixels. The total col- 
lection efficiency for the depletion region remains unchanged since the sum of 
the collection efficiencies in the overlap region is unity. Because of this over- 
lap, the MTF due to the collection process is reduced from Rp to Rp, where 


t» . wb . w( a + b) /t->\ 

Rp = smc-g- smc — 2 (12) 

V. CHARGE TRANSFER INEFFICIENCY EFFECTS 
ON A LARGE -AREA ARRAY 

In an image sensor such as a vidicon, the electron beam provides a direct 
readout of each pixel. On the other hand, in a CCD, the signal charge from 
each pixel must first be transferred to the output for readout. In practical 
devices, the transfer process is less than perfect, with some charge redistribu- 
tion and loss occurring in the transfers and the dispersion increasing with the 
number of transfers. The key characteristic of a CCD is its charge transfer 
efficiency. The frequency- dependent distortion introduced by charge transfer 
inefficiency for a CCD delay line for the case where the fixed loss is negligible 
has been calculated utilizing both Z transform (Refs. 5 and 6) and Fourier 
analysis (Refs. 7, 8, and 9). The transfer function, derived from the transform 
of the impulse response in the time domain, for such a delay line is given by 
(Refs. 5 and 7) 


where 


H(w) = A(w) exp [-i</?(w)J 


M ■ [— 


u - o 2 ; 

- 2e cos ud 


n/2 


(13) 
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and 


</?(w) = n tan 


e sin odd 
1 - « cos ud 


where e is the fraction of charge left behind per transfer, n is the number of 
transfers, and d is the pitch of the pixels. The modulus A(w) is the attenuation 
factor, and the argument is the phase shift. For « << 1, these expressions 

can be approximated by (Refs. 7 and 10) 


|a(w)| =s exp (-ne)(l - cos ud) 
</?(<*>) s nc(wd - sin wd) 


(14) 


In the time domain, which is the more suitable domain to calculate pulse 
train dispersion effects, the charge Qj^ntg) in a charge packet at a time ntg, 
transferred through N pixels, is given by (Ref. 11) 


Q N (nt 0 ) = (1 - e) N Q 0 [(n - N)t Q ] + Q R (nt Q ) 


(15) 


where 

^^CC) 

°R (nt 0 ) = (1 ' ° N Xj( N m m ) €mQ ot (n - N - m)t o3 

m= 1 ' • 


Here tg is the clock period, Qq^Iq) is the size of the charge packet inserted at 
the input at time ntQ, and is the contribution from the trailers from the pre- 
ceding charge packets. 

Utilizing Equations (14) and (15), the effects of charge transfer inefficiency 
on the MTF and point-spread function for a large-area CCD imager can be calcu- 
lated. The results are shown in Figures 14 and 15. Figure 14 is the calculated 
MTF for a 400 x 400 three-phase CCD imager with a 0.9 x 0.9 mil pixel size 
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currently under development (Ref. 1) for three locations in the image plane. 
Overlapping pixels with trapezoidal collection profiles are assumed as shown. 
Charge transfer efficiencies in the anticipated range for a buried- channel 
device were used. For f « *nyQ’ output modulation for a sine- wave input 
will be identical to the MTF, For optical inputs near however, phasing 

effects become significant and the output modulation will depend upon the phas- 
ing (Ax) between the optical input and the CCD, as shown in Figure 16, where 
the relations between the MTF and device response for sine- and square-wave 
inputs have been summarized (Ref. 12). 

Figure 15 is the corresponding charge dispersion at the same locations 
for various optical inputs. The first column shows the effects of charge trans- 
fer inefficiency on a single 1, which would roughly be the point-spread function 
for this array. For a single 1, the ratio of the first trailer to the primary 
charge packet is equal to ne. Column 2 indicates the output modulation at the 
Nyquist limit for the in-phase condition. Figures 14 and 15 indicate how the 
MTF and point- spread function can be expected to vary continuously across the 
image plane and represents the ideal case, where the charge transfer loss is 
uniformly distributed. In practical large-area devices, however, charge 
dispersion will usually occur because of the presence of isolated regions of 
poor charge transfer efficiency. For such a case, the results shown in Figures 
14 and 15 can still be applied by determining an effective ne product for the 
region in question. 

VI. CONCLUSIONS 

The MTF and point-spread function for a large-area CCD imager have 
been calculated. Both the theoretical model and the experimental results show 
that the MTF loss due to lateral diffusion in the frequency range of interest 
(f < ^yq) is minimal for a back-illuminated CCD thinned to 10 pm, with the 
size of the pixel and its geometric collection profile being the dominant factor. 
These results suggest that the MTF could be increased by a further reduction 
in device geometry before diffusion effects become significant. The key role 
of charge transfer efficiency has been shown. A charge transfer efficiency q 
of 0.9999/transfer or better will be necessary to minimize charge dispersion 
effects and obtain a device where the MTF and point- spread function are limited 
only by pixel size. 
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Figure 1. 




OPTICAL INPUT 


Figure 2. Thinned three-phase anodized AICCD (cross section along 
vertical channel) 
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Figure 4. MTF for lateral diffusion (L,p = 15 jim) 
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Figure 8. 
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Figure 9. Test setup for sensitivity profile measurements 
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RELATIVE OUTPUT 
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Figure 11. Sensitivity profile for three adjacent pixels (blue filter) (scan i. 
to vertical channel) 
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RELATIVE OUTPUT 
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Figure 12. Sensitivity profile for three adjacent pixels (blue filter) (scan II 
to vertical channel) 
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Rj-sinc (u)bl2) sine (oi(a + b)/2) 

Figure 13. Schematic illustrating overlapping depletion collection profile 
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Figure 14. Calculated MTF for a 400 x 400 three-phase CCD imager 
(MTF = n / n x MTF (nt) X R ) 
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Figure 15. Charge dispersion due to transfer inefficiency (linear model 
n + e = 1) 
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